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CN . Abstract 

>• ! We consider the time evolution of a discrete state embedded in a continuum. 

ly-^ ! Results from scattering theory can be utilized to solve the initial value problem 

O I and discuss the system as a model of wave packet preparation. Extensive use is 

O '• made of the analytic properties of the propagators, and simple model systems 

are evaluated to illustrate the argument. We verify the exponential appearence 
of the continuum state and its propagation as a localized wave packet. 
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This article will appear in a Special Issue of Journal of Modern Optics in 1997 

1 Introduction 

Quantum mechanics knows many examples of discrete states coupled to a continuum. 
Much work has been devoted to these problems in the past. In this paper we are going 
to review these works briefly, but reconsider them as methods to prepare time dependent 
states. All tunnelling problems are of this type, but in addition we have similar phenomena 
occurring when a photon excitation is followed by leakage into continuous outgoing states. 
A simple physical example is provided by autoionization, where an excited state is 
coupled to atomic ionization according to 

A + hu^ A* ^ A+ + €'. (1.1) 

In molecules, the analogous process is the predissociation reaction 

AB + huj ^ AB* ^ A + B, (1.2) 



and in chemical reactions the transition state may give a similar behaviour 

AB + C + hu}^ABC*^A + BC. (1.3) 

In a nucleus, ABC* would be a compound nucleus or a doorway state. In particle physics, 
the particle production processes provide additional examples. 

When the final states form a continuous spectrum, the decay of the initially prepared 
state is often described by an exponential. This was first derived from quantum theory 
by Landau [|T|, but the procedure is usually referred to by the names Weisskopf and 
Wigner |^. These features were found to appear exactly in the treatment of the Lee 
model 0, where, however, the primary interest was devoted to the renormalizability of 
the theory. The appearance of decay as a consequence of the analytic properties of the 
propagators became a central issue for the theoretical discussion after Peierls had pointed 
out that the functions needed to be continued to the second Riemann sheet 0]. These 
discussions can be found in Refs. and , and the results are summarized in the context 
of scattering theory in the text book . 

When several singularities are close to each other or multiple poles occur, there appears 
non-exponential time dependence as discussed by Mower pj; the corresponding problem 
for a generalized Lee-model is treated in Ref . |p . In molecules the high density of states 
gives rise to a similar physical situation through nonradiative transfer of excitation; for 



a review of this theory see Ref. |T0|. Many more works from this early period could be 
cited, but these may suffice for the present. 

The works above are mainly formulated as scattering problems and, except for the 
evaluation of probability decay, they tend not to look into the detailed time evolution of 
the processes. The use of well controlled laser pulses in the femtosecond time range has 



recently made it possible to follow microscopic processes in space and time |TT|. Thus 
one may resolve the Schrodinger propagation following a fast preparation of an initial 
state. It is then interesting to reconsider the theories above from the point of view of 
state preparation. If we excite a resolved initial state coupled to a continuum, which kind 
of wave packet can we prepare into the continuum? How is this state emerging and how 
is it propagating? 



This can be seen as a continuation of our earlier work on molecular dynamics [|T2 



and on electrons in semiconductor structures |T^. Here we consider explicitly the time 
dependence following an initial preparation. Our state functions can thus stay normalized 
at all times in contrast to the situation in a scattering description. Thus we can avoid 
paradoxes of exponentially growing amplitudes such as discussed by Peierls in Ref. [Q . 
The preparation of well localized wave packets is an interesting and challenging prob- 
lem in atomic physics. Various excitation processes have produced electronic wave packets 



on bound Rydberg manifolds; see e.g. Refs. JTSf and |T6|, but it is not clear if these can 



be launched into freely propagating states. Molecular dissociation, following coherent ex- 



citation from the ground state |jT2|, is expected to prepare localized wave packets in the 
reaction coordinate. In simple cases, this can be seen as propagation in the laboratory 
too. To prepare electronic wave packets in semiconductors poses huge experimental dif- 



ficulties [T7|; to retain their coherence seems almost impossible. Excitonic wave packets 
dephase more slowly, but their theoretical description is even more difficult. 

In this paper we adopt the presumption that an experimental procedure exists to 
prepare a pure isolated state. The interaction with the continuum is set in, and we can 
follow both the growth of the continuum states and the decay of the initial one. This 
makes it possible to consider the process as a state preparation. We can describe the 
disappearance of the initial state and the growth and propagation of the state prepared 
in the continuum. 

In order to achieve our goal, we postulate the validity of a simple model Hamiltonian. 
It is of the Lee model type and can hence be solved exactly to give formal expressions for 
all quantities of interest. In order to treat the initial value problem we utilize the Laplace 
transform, but in order to retain agreement with the conventional Fourier transform 
treatments, we denote the Laplace transform variable by —iu; the resulting transform is 
named the Laplace- Fourier transform (£jF-transform). 

We present the model in Sec. 2 and its formal solution in Sec. 3, where also the general 
behaviour of the solution is discussed. In Sec. 4 we present some simple models which 
illustrate the features of the general discussion. 

In Sec. 4.1 we assume that the continuum spectral density can be described by a 
simple complex pole, the generalization to several poles is straightforward. This approx- 



imation has acquired renewed interest after Garraway |]I8| has showed that the multiple 
pole approximation can be put into a Lindblad form amenable to a Monte Carlo simu- 
lation in terms of hypothetical pseudomodes, which describe the non-Markovian effects. 
However, he finds that consistency may require that the pseudomodes are coupled in their 
Hamiltonian. 

Section 4.2 presents a simple model where the role of analytic continuation can be 
elucidated and the exact behaviour of the wave packet in the continuum can be calculated. 
The role of the process as a preparation is clearly seen in this model. In Sec. 4.3 we 
try to construct a physical system having the properties of our earlier models. This 
is selected from the field of molecular physics ^2 and tests how weak coupling to a 



molecular continuum can be utilized to prepare a wave packet. We find that this is, 
indeed, possible, but the system has got many complementary features which we will 
discuss in a forthcoming paper. 



Section 5 presents a brief discussion of our results and the conclusions. 

2 The model calculation 

We are considering a model with one single state embedded in a continuum. This is 
described by the Hamiltonian 

H = Ho + V 

Ho = uo |0)(0| +Jde\e)e{e\ (2.1) 

V = Jde (V;|e)(0|+V;*|0)(e|) 

in an obvious notation. Observe that this is the form of a general tunnelling Hamiltonian. 
The initial state is supposed to be the isolated state 

|V^(t = 0))=|0). (2.2) 

Because we are going to consider time evolution, we introduce the Laplace-Fourier trans- 
form in the form 



oo 



^(tu) = -i dte"^'^lj{t) = CT{^p). (2.3) 



This gives us analytic functions in the conventional half-plane Im(c<;) > 0. We have the 
usual properties 

CJ^i^) = -i {coij{uj) - ij{t = 0)) 

jirn {ujij{uj)) = ip{t = 0) = ^^ (2.4) 

lim (utpiuj)) = ifjit = oo). 
The inversion is achieved by 

+oo+ia 

m = ^ I e-^-'^{u)dt, (2.5) 

— oo+ia 

where a > 0. With these definitions the solution of the Schrodinger equations is given by 

i>{u) = Giuj)^o, (2.6) 



where the propagator is 



G(.) = -^. (2.7) 



It is easy to verify the consistency of these relations. 

We are now going to partition the problem above in order to separate the time evolu- 
tion of the single state and that of the continuum. Thus we introduce the projectors 

P = |0)(0| 

(2.8) 
Q = ^-P = Ide\e){e\. 

Both projectors commute with the Hamiltonian Hq. 

The results that we need were derived long ago in scattering theory, but for easy 
reference we summarize them in the Appendix. We introduce the definition 

O^^ = XOY, (2.9) 

where O is any operator and X,Y = P or Q; the operators O^^ are denoted by O^ 
simply. 

We introduce a connected operator F which has the single state pole removed from 
the intermediate states by writing 

T = V + VG'^T, (2.10) 

where the unperturbed propagator is 

Goiuj) = ^— . (2.11) 

uj- Ho 

In terms of this operator, we can now write the exact solution for all the propagators 
needed in the form 

G^ = ((Go^)-i_r^)"' (2.12) 

G-Q^ = G^T/Q^G^ (2.13) 

G« = Go^ + G^(r« + r«^G^r^«)G^. (2.14) 



Thus we can obtain all desired partitioned propagators if we can solve the integral equation 
( [2.101 ). In the present case, this can be solved trivially if we notice that 

yP ^yQ = 0. (2.15) 



(2.16) 



From this it follows that 

YQp = yQP 

yp = i^^qg;?^«^_ 

This gives the solutions 

(0|G|0) = ^-- (2.17) 



The self-energy function is given by 



S(cu) = (0|r|0) = fde(^^] 



(2.20) 



where we have introduced the spectral density 

D{e) =1 V, |^> 0. (2.21) 

This notation is preferable because we can then include a density-of-states function in the 
spectral density D{uj) when needed. 

With the initial condition {^^ ) we find the amplitude of the isolated state to be 

-iHt I n\ ^ f „-i<^t 



Ait) ^ (0 I e--' I 0) = - / e— ' — - dt. (2.22) 

2tt J \uj — uJo — h[iO)l 

—oo+ia 

Because the integrated function is analytic in the upper half plane, the time evolution is 
determined by the singularities of the integrand on the real axis and below. We are thus 
looking for solutions of the equation 

h{uj) = uj-uJo-i:{uj) = 0. (2.23) 

As can be seen from Eq. ( ^.2(j| ) the function Ti{uj) has got a branch cut along the real 
axis. Its weight is given by D{uj). 



3 General properties of the solution 

From the definition (|2.20|) it follows directly that 

lmJ:{uj±iri) = TT^D{uj); (3.1) 

it is assumed that r/ ^ 0. Consequently, the function has a branch cut at every value of 
u such that D{uj) differs from zero. 

All physical systems have energies bounded from below, thus all integrals over the 
variable e must start at a finite value /i. If the energy of the isolated state lies below this 
spectral cut-off 

uo < /i, (3.2) 

the continuum only serves to renormalize the state. Then we may write 

E(.)=/* £^ ^=/*^_/*^(l) (,-..) + ,.. (3.3) 

This gives the propagator (|2.17] ) the form 

(0 I G I 0) = ^. , (3.4) 



where the state renormalization constant is given by 



Die) 
[ujQ - e) 



-1 



Z={l + jde^^^^^\ (3.5) 



and the renormalized energy is 

Coq = ulIq + Z de < uq. (3.6) 

J [ujo - e) 

These results are valid to order ( "°J^° j , and we still have ujq < jjl. The only effect of the 
continuum is to push the isolated level away and decrease the overlap between the initial 
bare state and the renormalized state to Z < 1. This state does not decay but continues 
oscillating with the frequency uJo forever. The reminder of the initial state is lost into the 
continuum as the coupling is switched on. 

When the isolated state is embedded in the continuum, uq > /i, it can decay. In the 
equation ( p.22| ) we require the zeros of the function h{uj) to be in the lower half plane. 
Let us see if such singularities exist by setting 

uj = uj' — ioj" (3.7) 



into h{uj). We find 

uj' - too" -ujo- fde ^ = 0. (3.8) 

Separating the real and imaginary parts of tliis equation sliould fix the oscillational fre- 
quency uj' and the damping uj". We write down the imaginary part as 

This equation obviously lacks solutions u" ^ and does not correspond to a physical 
result. This is obtained when we remember that the physical values of u approach the 
real axis from the upper half plane, and we should thus look for zeros when the function 
E(a;) is continued analytically into the lower half plane. This procedure will provide 
the singularities giving a contribution to the integral ( |2.22|) . We consequently need to 



continue T,{uj) across the real axis as discussed by Peierls ||T8[ 



The simplest manner to do this is to push the integration contour in Eq. ( p.22|) down 
below the real axis. Any isolated pole encountered in this manner is included as a single 
contribution to the time evolution in the manner shown in Fig. 1. For one single pole 
this works well, and the result is simple exponential decay as in the Weisskopf-Wigner 
approach. However, with several singularities only well isolated poles can be treated this 
way. When the poles are close to each other, this method cannot be justified, because 
its validity is based on a series expansion like that in Eq. (|3.3|) around each singularity; 
the radius of convergence can only be extended to the nearest singularity, and hence close 
lying poles interfere §]. Another limitation comes from the lower cut-off in the spectral 
weight D{uj). Because D{uj) = for cj < /i, this end point must be a singularity of the 
function. Thus there is a branch cut starting here, which can be moved around but must 
be pinned down at cu = /x. Poles too close to this cut-off cannot have a large radius of 
convergence and hence the branch cut cannot be neglected in their treatment. In particle 
theory these cuts derive from particle creation thresholds. 

In the region where D{uj) is an analytic function of a;, we can achieve the analytic 
continuation of the self-energy by defining the new function 

S+(cu) = j de ( ^^) - 27riD{uj). (3.10) 

Using (|3.1|) one can easily verify that the function S+(u;) has no singularity when the 
real axis is crossed from above. When this process works, it is a simple way to look for 
singularities on the second sheet of the lower half plane. 

When the result ( |3.1(J| ) is inserted into the equation ( |2.23| ) we find the relation 

' /■ , / X / uj' — e \ dD(uj') „ 

._,„_/,,fl(,)^__,__j+2,_i_J. =0, (3.11) 




(X 



\ 



X 



Figure 1: In the complex w-plane, the Laplace- Fourier transform is inverted by integration 
along the line [—00 + za, +00 + ia\. The analytic continuation can be effected by moving 
the contour into the lower half-plane. When a pole is encountered, the contour must circle 
it, but can continue down leaving the pole encircled by an isolated contour. 



where we have anticipated that the imaginary part is small. In addition we have 




u = 2mDiu; ) - / deDie) { — ^^;^ . (3.12) 






The values of the oscillational frequency lo and the damping lo are determined by the 
coupled equations ( |3.11| ) and ( p.l2| ). The oscillational part determined from ( |3.11| ) is 



discussed in detail by Cohen- Tannoudji ||2^. In the limit a; — > the results simplify. 
When the damping is small, we obtain 

J ^TxDiuj), (3.13) 

which directly provides a consistency check. The Weisskopf-Wigner result thus emerges 
in the weak coupling limit from our prescription for analytic continuation. 

It is instructive to realize that in the weak coupling limit, the prescription we have 
used for analytic continuation is to replace S(a;) just below the real axis by its value just 
above; i.e. we are performing an analytic continuation using only the zeroth order term 
in a Taylor expansion across the real axis. Thus the analytically continued equation ( p.8|) 
becomes 

J - iuj" -uq- I de — = 0. (3.14) 

For small values of u" this gives back all results derived previously. It is to be noted that, 
inside the integral, the precise value of uj" does not affect the result, and we can obtain 
the correct damping using the conventional infinitesimal irj prescription. Our discussion 
is in fact a derivation of this rule. 

An instructive way to consider the equation ( |3.14| ) is to regard the imaginary part in 



the denominator as a small dissipative part of the energies in the continuum 

e — ioj" ^ e — ie . (3.15) 

This provides an interesting physical interpretation of the results. Adding any minute 
dissipative mechanism to the continuum we will find that the isolated state will decay with 
the rate (|3.13|) independently of the details of the dissipation of the continuum. Thus any 



system coupled perturbatively to a dissipative probability sink through a continuum will 
decay with a rate determined by the coupling strength and not by the actual dissipation of 
the continuum. In the perturbative limit, the probability flow into the continuum will be 
the bottleneck. The dissipated continuum acts as a proper reservoir unable to return the 
probability once it has received it. This interpretation of the procedure directly provides 
the right sign in the denominator of ( 3.14| ) to effect the proper analytic continuation. 
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When we calculate the integral in Eq. ( p.22| ) , we distort the contour into the lower half 
plane until we encounter the pole defined by (|3.11| ) and separate a contour cq circling this. 
If no other singularity were encountered, we could pull the contour to negative imaginary 
infinity and find no corrections. In a physical system this is, however, impossible because 
the spectral density D{uj) has necessarily got the lower limit fi as explained above. The 
branch cut starting here is pulled down to —ioo in the manner shown in Fig. 2. This 
contributes an additional contour ci as shown in the figure. Thus we find 



A{t) 



2^ 



dte 



-iujt 



1 



ci 



h (t) + exp 



-lUJ t — UJ t 



SM, 



dte 



-iuit 



UJ ^ iuj' 



LU 



(3.16) 



The second term is as expected in a perturbative calculation of the Weisskopf-Wigner type, 
and we have taken the residue Z of ( p.5|) to be unity. The decay rate of the probability 
is given by the well known expression 



'y = 2uj = 27vD{uj). 
However, the corrections from the cut are given as 



hit) 



27r 



<t 



[UJ 



ujo — Tij) (oj — ojq — T, 



d^, 



II 



(3.17) 



(3.18) 



where the infinitesimal parameter r] ensures that the self energy is evaluated on the two 
sides of the cut. From this form we can see that in the long time limit t — *> oo, only the 
value ^ ~ contributes. This would give no contribution without the singularity of E 
across the cut. Thus we can set ^ = everywhere except in the integral 



/•oo 

S/(yU - i^ + ?]) - S//(yU - i^ - ?7) = / de D{t 

J a 



1 



1 



/i 



ii + T] 

1 



/i 



1 



i + P + ITj ^ + u — if] 



= diy D{^ + iu) 



= -2mD{fi + i^). 
Because lj = fi is a threshold, we expect that we have 

for small values of ^. Inserting this into the equation ( |3.18| ), we find the leading term 

1 



(3.19) 



(3.20) 



h(t) = — ^Fa + l, , ,, 



(3.21) 
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Figure 2: Because physical spectra have a lower limit /i, there must be a branch cut 
starting at this value. When isolated poles are encircled by closed contours, here Cq, the 
branch cut can be drawn from /i to — cxd, thus contributing another contour Ci. 

This type of asymptotic result was given for scattering theory already in Ref. 0. When 
the exponential in ( 3.16 ) has decayed, the power dependence ( p. 21 ) will dominate for long 
times. 

The treatment given above does, of course, assume that there are no additional poles 
encountered. If these are well separated, they contribute simply additional exponential 
terms, but for poles situated close to each other interferences occur; for a discussion of 
these cf. the discussion by Mower M. 



4 Simple models 



4.1 A pole in the spectral density 

Many systems display a moderately sharp spectral feature in its density of states. Hence 
one may often approximate this with a Lorentzian shape, which simplifies the treatment 
considerably. Such assumptions are often utihzed to describe scattering from condensed 
matter [^], but it is useful in theoretical discussions too. Recently Garraway has utilized 
it to describe non-Markovian effects in terms of pseudomodes |T^ . 
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We assume that the spectral function has the simple shape 

DU) ^ 



it - aY + 62 



(4.22) 



This satisfies the positivity requirement, but it lacks the lower limit necessary in physical 
systems. Thus no power law terms are expected. 

With the expression ( |4.22| ), the analytic continuation to the second sheet is trivial. 
Inserted into Eqs. (|2.2CI|) and ( p.l7|) it gives the results 



S(^) 

(0 I G{uj) I 0) 



nA 



uj — a + ib 

{uj - a + ih){LO - ujq) - ^ 



(4.23) 



This form is found to satisfy the correct initial condition 

lim cj(0 I Giuj) I 0)^(0) = ^(0). 

w— >oo 

The general solution of the time evolution becomes 



A{t) = (i?i e*^+* + i?2 e 



in^t 



(4.24) 



(4.25) 



where Q± are the two roots of the denominator in Eq. ( [4.23|) ; from the form of the equation 
it follows that the sum of the residues is unity: Ri + R2 = I- 

Inserting the ansatz (|3.7|) into the denominator of Eq. ( 4.23 ) we find that the imaginary 
part follows as the solution of the equation 



'nA^ 



u 



UJ 



UJ 



ih 



UJ 



ii\2- 



(4.26) 



Because the right hand side of this equation is positive for tu < h there must be a 
solution for < cu < 6, irrespective of the value of u . This proves that both roots Vt± 
have negative imaginary parts and the solution A{t) decays to zero as we expect. Thus 
the pole approximation is consistent and only fails to reproduce the long time behaviour 
deriving from the threshold in the spectral density. In a similar manner the appearance 
of several poles can be discussed. Of special interest is the case when a pole with negative 
weight is encountered. This case is discussed in detail by Garraway [T^ . 
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4.2 The branch cut model 



In order to illustrate the concepts introduced in Sec. 3 we introduce a model where the 
Weisskopf-Wigner result emerges as an exact consequence. It is then possible to follow 
how the analytic continuation to the physical sheet works in detail, and the correction 
terms can be evaluated explicitly. 

The model is based on the following spectral density 

( A^; \e\<L 
D{e) = (4.27) 

[ 0; |e|>L. 

It is taken to be an essential feature of the model that L is taken to be large (infinite) in 
the final expressions. For finite L, it is trivial to compute the self-energy from ( |2.20D to 
be 



S(cu) = A^ log 



-{L + cv) 



(4.28) 



In this model we can explicitly see that the sign of the imaginary part clearly depends on 
the choice of branch of the function involved. We use the discussion in Sec. 3 to introduce 
the analytic continuation ( |3.1CI| ) 



If we now introduce the ansatz 



A' 



we find that 



where we have 



log( ;^J +Hf- -v+), 



and 



R. 



^T 



1 + 



„\2 



log ( -y^ ^ \ -2ni . (4.29) 

to = u — iuj . (4.30) 

log l ~i^^..x I = log i^] +^{^-- V^) , (4.31) 

(4.32) 

arctan I - — - — — 1 . (4.33) 

{L±u})J 

These relations are illustrated in Fig. 3. It is now easy to see that independently of 
whether we choose the cut of the logarithm function from the origin to +cxd or to — cxo, 
the limit becomes 

lim S+(c.;) = -inA^, (4.34) 

L— »oo 



Hf) 



\i^'i) +{i)' 



UJ 
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giving exactly the decay constant 



7 



27rA^ 



(4.35) 



cf. Eq. ( p.l7|) . The frequency shift vanishes in this hmit. Our treatment is not mathe- 
matically rigorous, and the contributions from the branch cut are not evaluated. They 
disappear when the lower spectral limit —L goes to infinity. It is, however, possible to 
evaluate the correction terms to order L~^ and verify the over-all validity of the picture we 
have developed. The model can also be extended to treat an asymmetric spectral weight 
function defined on the interval [L_,L_|_]. Here the limit appears slightly differently and 
modifications of the results above can be found. 







© 


R_e'^" 


cj"/L 


R+e'^ 


^1 




1 
1 



Figure 3: The branch-cut model gives the self energy S in terms of two vectors 
R± exp («v^±) in the complex plane as shown. 



With the result ( |4.35| ) we find an exact exponential disappearances of the amplitude 
of the initial state in Eq. ( [^.221 ). It is, however, instructive to evaluate the state in the 
continuum part too. The probability of emergence of this state cannot grow faster than 
the isolated state disappears, and consequently we expect it to need a time of the order 
7~^ to appear. From Eq. ( p.l^ ) we find for the continuum wave function 



^ 



cont 



>,t) 



de<t),{x){e\ G{uj) \ 0) e 



'iujt 



idui 
~2^ 



15 



$(x,t) 



K 



X) 



e -cuo + ^2 



) (»-*■' - 



g-ia;otg-7«/2^ 



(4.36) 



-ote-^*/2$(2,^0). 



Here 0e(x) is the continuum eigenfunction corresponding to the eigenvalue e, and $(x, t) 
is the state emerging from the process on the corresponding Hilbert space. This has got 

the spectral distribution 

7/27r 



Pi 



(4.37) 



as known from the Weisskopf-Wigner calculation. In addition, we can see how the con- 
tinuum state grows from an initial zero value to the emerging wave state $(x, t), which 
will travel according to the dynamics in the continuum as determined by the spectral 
distribution ( ^.37] ). This feature has usually not been discussed in the ordinary decay 
problems treated in the literature. 

In the present model the analytic behaviour can be seen directly. This is useful if 
we want to gain understanding of the features playing the central role in the general 
treatment. We can also calculate the emergence of the outgoing wave packet explicitly 



([4.36|) . In order to see the shaping of this wave packet we consider one more model where 
the features discussed can be followed numerically. 



4.3 A wave packet model 

As a model we choose a coupled pair of energy levels with one harmonic potential and 
one potential slope providing a continuum. The Schrodinger equation is 



.d_ 

''di 



^1 

1p2 






1„2 
2 



V 



y (-fe-/5^+7i 

This corresponds to a scaled molecular problem with the parameters 

V2. 



i^2 



m 



Because the zero point energy is given by 

1 



-ui 



UJ 



1 

71' 



(4.38) 



(4.39) 



(4.40) 



the ground state energy of the harmonic oscillator coincides with the continuum energy 
at the origin; see Fig. 4. To the extent we can neglect the influence of the first excited 
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Figure 4: The molecular energy levels have one harmonic potential and one straight slope. 
The parameters are chosen such that the energies of the two levels coincide at x = 0. 

harmonic oscillator state, we can consider this model as a realization of a situation where 
the initial oscillator ground state 

„2 \ 



M^,t = 0) = {nV2)-'/^ 



X 



exp 



2V2, 



(4.41) 



is embedded in the continuous spectrum of the slope; the coupling is given by V. When 
this is in the perturbative regime 

— « u; (4.42) 

we expect the Weisskopf-Wigner treatment to hold and the initial state to decay nearly 
exponentially. Setting 

V = 0.5; f3 = 3.0 (4.43) 

we satisfy ( [4. 421 ) and obtain the decay of the initial state 



Piit) 



ipi{x,t) p dx 



(4.44) 



shown in Fig. 5. As seen, there is a clear range of times over which the initial state decays 
exponentially as expected. By looking at the state emerging on the second energy level 
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ip2{x,t) shown in Fig. 6 we can see that the emerging state hke (|4.36|) indeed produces 
a wave packet which travels down the potential slope while it spreads according to the 
requirements of quantum mechanics. In the present model, however, the initial state 
is not depleted totally by the exponential decay. This derives from the fact that the 
crossing potential surfaces form quasibound states in the adiabatic representation, and 
the population trapped in these states can then ooze out only over times much longer 
than those characterizing the exponential decay. 

The model calculation presented in this section shows that our system ( |4.38| ) can serve 
as an illustration of Weisskopf-Wigner decay, and it proves that such decay can serve as 
a wave packet preparation device. However, we have found that even this simple system 
works only when the parameters are chosen right, and it contains features not expected 
from the simpler models discussed above. We consider these aspects in a forthcoming 
publication, where we shall present the general behaviour of the harmonic state resonantly 
coupled to a sloped potential. Here we have only presented a single case as an illustration 
of our general theoretical considerations. 



5 Discussion 
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Figure 5: The system is prepared on the ground state of the harmonic potential, and the 
probability that it remains there is plotted as a function of time. The probability leaves 
the level in an exponential manner as this semilogarithmic plot shows. 



In this paper we have reviewed the scattering theory approach to partitioning of the 
Hamiltonian time evolution in a quantum system. We are especially interested in the 
situation where one (or more) discrete states leak into a continuum, i.e., decay of qua- 
sistationary states. By choosing a Laplace-Fourier transform in time, we concentrate on 
solving an initial value problem and following the system as it irreversibly transfers its 
probability into the continuum. This is regarded as a wave packet preparation procedure. 

The emphasis on time evolution has been motivated by the recent experimental progress 
in pulsed laser technology. It is now possible to excite a single state selectively, control its 
coupling to other states, including continuum ones, and follow the ensuing time evolution 
of the quantum states. Such work has experienced tremendous progress in molecular in- 
vestigations recently, but also time resolved spectroscopy of semiconductor structures is 
possible. 

Starting from a giving initial time, the future evolution of quantum systems is deter- 
mined by the analytic properties of the propagator functions in the complex frequency 
plane. In order to obtain the correct physical behaviour these have to be continued an- 
alytically to a second Riemann sheet in the lower half plane. Any pole encountered will 
give rise to a resonance but the physical behaviour must also receive contributions from 
unavoidable branch cuts. 

The necessary analytic continuation introduces simple relations between the self- 
energy function in the upper half-plane and that in the lower. In the weak coupling 
limit, the Weisskopf-Wigner case, the prescription for continuation is found to be equiv- 
alent with adding a small dissipative mechanism to the states in the continuum. The 
causality requirement on the sign of this dissipation automatically affects the necessary 
analytic continuation. Physically we can understand this as follows: When probability 
leaks into the continuum through the weak coupling link, it is dissipated by all modes 
and oozes into the unobserved degrees of freedom providing the continuum damping. In 
scattering theory this mechanism is simply the outgoing boundary condition at infinity; 
all scattered modes are irreversibly lost at the edges of our universe. The bottle neck is 
provided by the weak link, and the rate of loss of the initial state is fixed by this. The 
dissipation in the continuum is only present to prevent any return of probability. Hence 
its details are inessential to the decay, its presence however is necessary. 

We have discussed the analytic behaviour of the propagators in detail, and illustrate 
the properties by simple model calculations. Most of the treatment follows the tradition 
in this field and carries out the argument in the energy representation. However, to 
prove the existence of the phenomena discussed, we display the behaviour of a model 
describing laser coupling of a bound molecular state into a dissociating continuum. In 
addition to showing the expected exponential decay, the model system proves that the 
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P 2 




Figure 6: The probability \ip2{x,t)\'^ leaving the bound oscillator state appears on the 
sloping continuum potential. Here we can see it appear and fall down the slope in the 
form of a spreading wave packet. Part of the probability gets captured near the position of 
the original bound state because of trapping in the potential well formed by the adiabatic 
levels. This contains a component on the continuous manifold, which explains the trapped 
portion. Even this must eventually leave the trap but clearly over a much longer time 
scale. 



state prepared on the continuum is in the form of an outgoing wave packet. The model, 
however, also displays features not simply describable in the model calculations performed 
in this paper; we will return to a detailed exploration of the decay characteristics in a 
forthcoming publication. 

The purpose of this paper is to draw attention to the necessity to investigate the 
time evolution in coupled quantum systems. This provides a broad range of interesting 
quantum models, which may also shed light on the phenomena occurring in dynamical 
processes initiated and explored by pulsed lasers. 
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Appendix 

The algebraic manipulations we need in this paper have been provided by scattering 
theory long ago 0. Here we summarize the main results only. 

The perturbation series for the propagator is generated by the relation 

G = Gq + GqVG. (A.l) 

Introducing the operator T by the relation 

VG = TGo, (A.2) 

we find the equation 

G = Go + GqTGo. (A.3) 

The operator T is obtained from the equation 

T = V + TGqV. (A.4) 

The disconnected operator F is defined by the relation 

T = V + VG^T, (A.5) 

where the projector Q = 1 — P has been introduced. The missing singularities in the 
intermediate states can be reintroduced to get the T operator 

T = r + TG^T. (A.6) 

When we remember the projectors commute with H, we can easily solve 



G^ = Gni + T^G^ 



and 



K)"'-r^, 



(A. 



In the same way we can calculate G"^^ from Eq. ([A. 1|) and G"^ from Eq. (|A.3|) . The results 
are given in Eqs. ( p. 13 ) and ( 2.14 ). 



These results are all exact, and they can be used to solve for the propagators of any 
Hamiltonian split into Hq and V and partitioned by projectors P and Q which separate 
the Hilbert space into subspaces not coupled by V . 



21 



References 

[1] L. D. Landau, Z. Physik. 45, 430 (1927). 

[2] V. Weisskopf and E. Wigner, Z. Physik. 63, 54 (1930). 

[3] V. Glaser and G. Kallen, Nuclear Physics, 706 (1956). 

[4] R. E. Peierls, p. 296 in Proceedings of the 1954 Glasgow Conference (Pergamon Press, New 
York, 1955). 

[5] B. Zumino, p. 27 in Lectures on Field Theory and the Many-Body Problem, ed. E. R. 
Caianeho (Academic Press, New York, 1961). 

[6] M. Levy, p. 47 in Lectures on Field Theory and the Many-Body Problem, ed. E. R. Caianello 
(Academic Press, New York, 1961). 

[7] M. L. Goldberger and K. M. Watson, Collision Theory, Chapter 8 (J. Wiley & Sons, New 
York, 1964). 

[8] L. Mower, Phys. Rev. 142, 799 (1966). 

[9] P. V. Ruuskanen, Nuclear Physics B22, 253 (1970). 

[10] J. Jortner, S. A. Rice and R. M. Hochstrasser, p. 149 in Advances in Photochemistry, Vol.7, 
eds. J. N. Pitts, Jr., G. S. Hammond, and W. A. Noyes, Jr. (Interscience, New York, 1969). 

[11] M. Gruebele and A. H. Zewail, Physics Today, p. 24, May 1990. 

[12] B. M. Garraway and K.-A. Suominen, Rep. Prog. Phys. 58, 365 (1995). 

[13] S. Stenholm and M. Kira, Acta Physica Slovaca 46, 325 (1996). 

[14] R. Peierls, More Surprises in Theoretical Physics, Sec. 5.3 (Princeton Univ. Press, Prince- 
ton, 1991). 

[15] G. Alber and P. Zoller, Phys. Rep. 199, 231 (1991). 

[16] M. Nauenberg, C. Stroud, and J. Yeazell, Scientific American, June 1994, p. 24. 

[17] C. Weissbuch and B. Winter, Quantum Semiconductor Structures, (Academic Press, New 
York, 1991) . 

[18] B. M. Garraway, to be published. 

[19] B. M. Garraway, private communication. 

[20] C. Cohen- Tannoudji, J. Dupont-Roc and G. Grynberg, Atom-Photon Interactions: Basic 
Principles and Applications, p. 239 (Wiley-Interscience, New York, 1992). 

[21] P. C. Martin, p. 37 in Many-Body Physics, Les Houches 1967, eds. B. DeWitt and R. 
Bailian (Gordon and Breach, New York, 1968). 



22 



